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Using the full set of the LHC Higgs data from the runs at 7 and 8 TeV center of mass energies 
that have been released by the ATLAS and CMS collaborations, we determine the couplings of the 
Higgs particle to fermions and gauge bosons as well as its parity or CP composition. We consider 
ratios of production cross sections times decay branching fractions in which the theoretical (and some 
experimental) uncertainties as well as as some ambiguities from new physics cancel out. A fit of both 
the signal strengths in the various search channels that have been conducted, H — > ZZ, WW, 77, rr 
and bb, and their ratios shows that the observed ~ 126 GeV particle has couplings to fermions and 
gauge bosons that are Standard Model-like already at the 68% confidence level (CL). From the 
signal strengths in which the theoretical uncertainty is taken to be a bias, the particle is shown 
to be at most 68% CP-odd at the 99%CL and the possibility that it is a pure pseudoscalar state 
is excluded at the 4<r level when including both the experimental and theoretical uncertainties. 
The signal strengths also measure the invisible Higgs decay width which, with the same type of 
uncertainty analysis, is shown to be r^/Pf^ < 0.52 at the 68%CL. 

1. Introduction 

The ATLAS and CMS collaborations have released their analyses of the mass and the pro- 
duction times decay rates of the 126 GeV Higgs-like particle using the full set of data collected 
so far, « 5 fb" 1 at v^s = 7 TeV and « 20 fb" 1 at ^ = 8 TeV [M]. This closes a very 
successful first run at the LHC, which culminated with the historical discovery of the state in 
July 2012 [51 16]. To be convinced that the observed particle is indeed the Higgs boson that 
is responsible of the spontaneous breaking of the electroweak symmetry [7J [8], one needs to 
prove that the particle: i) has spin-zero, ii) is a CP-even state, Hi) couples to fermions and 
gauge bosons proportionally to their masses and, ultimately, iv) has a self-coupling that is also 
proportional to its mass. While there is little doubt on the spin-zero nature of the observed 
state 1 , and the probing of the self-coupling has to await for a high-luminosity LHC [IT] or 
a future lepton collider, a first determination of the couplings |12hT^] and the CP-properties 
[151 116] can be performed with the current results. The ATLAS and CMS collaborations them- 
selves have given a first "portrait" of the observed particle which indicates that indeed it has 
the properties of a Higgs boson and even more, the properties of the unique Higgs boson that 
is predicted in the Standard Model (SM) [IH3|. 

However, in the case of the particle coupling determination, the ATLAS and CMS analyses 
suffer from two serious drawbacks [TTJ [JS]. The first one is that the signal strength modifiers 
Hxx-, that are identified with the Higgs cross section times decay branching ratio normalized 
to the SM expectation in a given H — > XX search channel, are affected by large theoretical 
uncertainties that are now becoming a dominant source of error. Indeed, the combined theo- 
retical uncertainty in the rate of the by far dominant Higgs production process at the LHC, 
gluon fusion gg — > H, is estimated to be of order ±15-20% [El GO] even before it is broken 
into jet categories which significantly increases the uncertainty [211 [22] . The uncertainty is 



[1] The observation of the H —> 77 decay rules out the spin-1 case [9] and the graviton-like spin— 2 possibility is extremely unlikely 
and, from the particle rates, is ruled out in large classes of models I10| . 
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similar in the vector boson fusion channel when the large contamination from the gg — > Hjj 
process is taken into account [22]. Another drawback of the analyses [T71 [18] is that they involve 
strong theoretical assumptions on the total Higgs width since some contributing decay channels 
not accessible at the LHC are assumed to be SM-like and possible invisible Higgs decays in 
scenarios beyond the SM are supposed not to occur. 

In this letter, we consider ratios of Higgs production cross sections times decay branching 
fractions in which these two sources of uncertainties and, eventually, also some systematical 
and parametrical uncertainties such as the error on the luminosity measurement and the one 
on the Higgs branching ratios 2 [201 123], should be absent [T71 [IS] . Using the Higgs signal 
strengths in which the theoretical uncertainty is taken to be a bias and not a nuisance (as is 
done generally), as well as the ratios /i 77 / fizz and fi TT j fiww which are free from the ambiguities 
above, we perform a fit of the latest ATLAS and CMS data and conclude that, already at the 
68% confidence level (CL), there is no deviation of the Higgs couplings to fermions and gauge 
bosons from the SM expectation. 

On an other front, the attempts made so far for the determination of the CP nature of the 
particle mainly exploit the kinematical features of the H — > VV decays with V = W, Z [15J 
or the production in the vector boson fusion V*V* — > H and Higgs-strahlung V* — > VH 
processes [IB] . Since a CP-odd particle has not have tree level couplings to VV states 3 , all 
these processes project out only the CP-even component of the HVV coupling [27] and the 
considered distributions can be thus only those of a ++ state. 

We will show that a much better way to measure the CP composition of the observed Higgs 
state is to consider the measured signal strength in the H — >• VV decays. Using fizz, we 
demonstrate that, if the magnitude of the Higgs couplings to fermions is as in the SM, the pure 
CP-odd possibility is excluded at the 4cr level, irrespective of the (mixed CP) Higgs couplings 
to light fermions. 

Finally, assuming that the Higgs couplings to fermions and gauge bosons are SM-like and 
using mainly the signal strength f^zz, one obtains a limit on the rate for invisible Higgs decays, 
r^ v /r| M < 0.52 at the 68%CL (see also Ref. pU [H] for an indirect limit), that is stronger 
than the one obtained from direct invisible Higgs searches [231 |2"B"] . 

2. The signal strengths and their ratios 

Let us first summarize the LHC Higgs data collected in the 2011 and 2012 runs for the various 
SM Higgs decay channels that have been searched for by the ATLAS and CMS collaborations: 
H ->■ ZZ* -> 4£± H -> WW* 2£2v, H 77, H ->• t+t~ and H ->• bb; we will 
ignore the additional search channels H — > fi + fi~ and H — > Z7 for which the sensitivity is 
still too low. In most cases, the various Higgs production channels have been used: the by 
far dominant gluon-gluon fusion mechanism gg — >• H (ggF) that has the large production 
rates but also the subleading channels, vector boson fusion (VBF) qq — > Hqq and Higgs- 
strahlung (HV) qq — > HV with V — W,Z] the top quark associated pp — > tiH mechanism 



[2] There are parametric uncertainties that affect the hadronic Higgs decay widths. For a f 26 GeV SM Higgs boson, this translates 
into an uncertainty on the total Higgs width which is found to be of order 4% in Ref. 1231 and slightly higher in Ref. 1201 . 

[3] The effective VV coupling of a pseudoscalar A boson, <x M^V^V^y with V^v = e M " pCT V r p (T , should be generated through tiny 
loop corrections. To make this coupling as large as the SM tree— level HVV coupling, one needs a very low new physics scale 
that would spoil the precision electroweak data. 
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(ttH) has too low a cross section to be relevant. At least the ggF and VBF channels have been 
considered in the H — > ZZ, WW, 77 and H — > t + t~ channels, while in the case H — > t + t~ 
and H — > WW* — > Ulv decays, also the HV production mode in which the H — )■ bb decay has 
been searched for, has been considered. 

ATLAS and CMS have provided the signal strengths for the various final states with a lu- 
minosity of, respectively, pa 5 fb _1 for the 2011 run at ^/s = 7 TeV and pa 20 fb _1 the 2012 
run at y/s = 8 TeV. We will identify these /z values with the Higgs cross section times decay 
branching fractions normalized to the SM expectation and, for the H — > XX decay, one would 
have indeed in the narrow width approximation, 

_ a(pp -» H -)■ XX) a(pp ^H)x BR(H -)> XX) 

™ th " a(pp ^H^ XX)\ SM ~ a(pp H)\ SM x BR(H -> XX)\ SM ' ' ' 

From the experimental point of view (and in our fits), this would correspond to 

Arevts 

™c X p " e x a(pp H)\ SM x BR(H ^XX)\ SM xC> { ) 

where N^x stands for the measured number of events in the H — > XX search channel, e 
denotes the selection efficiency and C is the luminosity. 
In this paper, we consider the decay ratios D X x discussed in Ref. [T7] and defined as 

a p (pp H XX) _ a p (pp H) x BR(H -»■ XX) _ T(H -»■ XX) 
xx ~ aP{pp ^ H ^ VV) ~ aP{pp H) x BR(if -> VV) ~ f(ff ->• VV) ( ^ 

for a specific production process p = ggF, VBF, VH or all (for inclusive production) and for a 
given decay channel H — > XX when the reference channel H — > VV with V = W or /and Z 
is used. In these ratios, the cross sections a(pp — » H) and hence, their significant theoretical 
uncertainties will cancel out as discussed previously, leaving out only the ratio of decay branch- 
ing fractions and hence of partial decay widths. Thus, the total decay width which includes 
contributions from channels not under control such as possible invisible Higgs decays, do not 
appear in the decay ratios D xx . Some common experimental systematical uncertainties such 
as the one from the luminosity measurement as well as the uncertainties in the Higgs decay 
branching ratios also cancel out. We are thus, in principle, left with mostly the statistical 
uncertainty and some systematical errors 4 . 

The ratios Dxx involve, up to kinematical factors, only the ratios |cx| 2 /|cy| 2 of the reduced 
couplings of the Higgs boson to the particles X and V compared to the SM expectation, 
cx = gnxx/gnxx- dxx- For the SM Higgs boson with a mass of Mh = 125 GeV, the 
kinematical factors can be straightforwardly obtained using the program HDECAY [30J for the 
evaluation of the Higgs branching ratios when the SM inputs recommended by the LHC Higgs 
working group [19] are adopted; they are given in Ref. [17] for the various normalisations (and, 
for simplicity, are then set to unity in that paper). 



[4] The theoretical and common systematical uncertainties will completely cancel out only when the same selection cuts are applied 
for the different final state topologies in a given production process (the selection efficiencies should be = key, k being a 
constant). This is obviously not the case in all the channels that we are considering here. We will assume, nevertheless, that 
this will be the case and we consider that the remaining uncertainty in the ratio is, to a good approximation, only of statistical 
nature. We hope that in the future, with the much larger data sample that is expected, the ATLAS and CMS collaborations will 
analyze the various search channels under the same experimental conditions. 
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In practice, to take into account the fact that there are four different Higgs production 
channels with different topologies, the cross section part is more involved and the ratio Dxx 
can be more precisely written as 

eg a(gg^H) + e x ^a(Hqq) + e x v a(HV) + efajttH) rg^L m 
XX e^a(gg^H) + e^a(Hqq) + e* v a{HV) + ef H a{ttH) r g^ M 1 ) 

where the e x , provided by the ATLAS and CMS collaborations, denote the experimental effi- 
ciencies to select the Higgs events in the gg, VBF, HV, tiH production and H — > XX decay 
channels (exclusive cut categories are also considered). Dxx is only proportional to the above 
expression due to the presence of another identical ratio of cross sections but within the SM 5 . 
Nevertheless for almost the same selection efficiencies, and even for e x = key with k a constant, 
the expression eq. Q simplifies to that in eq. ([5]). 
In fact, the decay ratios Dxx can be simply written in terms of the signal strengths 

n ^XX ^ a(pp^H)|sMxBR( J H-^XX)|sM = BR(g^XX)| SM = r(^!X)| SM _ |gx| , , 

n vv ~ <r(pp->H) x BR(g-> VV ) BR(H->VV) F(H-yVV) | c 12 W 

^ vv a(pp->H)\ SM xBR{H->VV)\ S M BR(i2->.VV)| SM T(H^VV)\ SM 1 V] 

where have used as normalisations in the ratios the channels 6 H — > VV with V = Z or W. 

Nevertheless, performing the ratios of signal strengths leads to a loss of information and, in 
some case, we will need at least one signal strength to set the normalisation. Rather than using 
the global /x tot value obtained by combining all search channels, we will consider the cleaner 
H — > ZZ channel alone as it is fully inclusive and thus does not involve the additional large 
scale uncertainties that occur when breaking the gg^-H cross section into jet categories 7 . The 
combination of the ATLAS and CMS data in the ZZ channel gives 

fizz = 1.10 ± 0.22 cxp ± 0.2 th (6) 

where the first uncertainty is experimental and the second one theoretical and that we assume to 
be, conservatively, A th = ±0.2. It has been advocated in Ref. [TTt 120] that, since the main effect 
of the theoretical uncertainty (which has no statistical ground) is to modify the normalisation 
of the SM cross section, it should be considered as a bias (rather than a nuisance as in the 
case of the experimental error) and, hence, one needs to perform two separate fits: one with 
Hzz + A th and another with [izz — A th . 

3. A combined fit of the Higgs couplings 

In order to study the Higgs at the LHC we define the (now usual) effective Lagrangian, 

C h = c w guww H W+W-" + c z g HZZ H Z ^ (7) 
- c t y t Hi L t R - c c y c Hc L c R - c b y b Hb L b R - c T y T Hf L r R + h.c. 



[5] Without derivative Higgs couplings, the kinematics and selection efficiencies are as in the SM. 

[6] In fact, one can assume custodial symmetry and use the combined H — ¥ WW and H — > ZZ channels as a reference to increase 
the statistical accuracy of the normalization factor. The ratio Dzz = Czz/fw has been measured for instance by the ATLAS 
collaboration to be Dzz — \ c z\ 2 /\ c w\ 2 — l-G^o's El (with the error expected to be only statistical), hence, supporting this 
approach. 

[7] In addition, contrary to the global signal strength fitot, it does not involve the channel T(H — > 77) which, at least in the ATLAS 
case, deviates from the SM prediction and might indicate the presence of new physics in the H-fy loop which is not reflected in 
the other couplings. 
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where yt, c ,b,r = m t,c,b,r/ v are the Yukawa coupling constants in the mass eigenbasis (L/R 
indicates the fermion chirality and we consider only the heavy fermions that have substantial 
couplings to the Higgs boson), gnww — ^^wl v anc ^ 9hzz = M^/v are the electroweak gauge 
boson couplings and v is the Higgs vacuum expectation value. The c parameters are all defined 
such that the limit c — > 1 corresponds to the SM case. For the present task, we assume no or 
negligible new contributions to the Higgs couplings to photons or gluons, e.g. as induced by 
new particles. 

We will present the results for the fits of the Higgs signal strengths, /i« (i labels each channel 
and cut category investigated), in the plane Cf versus cy. We have chosen universal coupling 
corrections, Cf = c t = c c = = c T and cy = cw = cz, for an illustrative purpose. All the 
Higgs production/decay channels are considered here and the data used are the latest ones and 
are borrowed from Ref. [3T] for the Tevatron, from Refs. (32H2S] and the combined results of 
Refs. [HE] for the ATLAS collaboration as well as Refs. [3"] IH I3TH1T] and the combined analysis 
of Ref. [3] for the CMS collaboration and finally the LHC results presented at the Moriond 
conference this month [261 [29] . 

We will closely follow the procedure of Ref. [13] (to which we refer for the relevant details) 
for the fit of these data, with one major difference though, the treatment of the theoretical 
uncertainty. The errors used in the present fit are those given by the collaborations as quoted 
above and, thus, contain the experimental and theoretical uncertainties added in quadrature, 



SfXi = y <5/^|e X p + (fy-tilth! which completely dilutes the effect of Sfii\ t h- In the present analysis, 

we treat the theoretical uncertainty (and not as if it were associated with a statistical 

distribution) and perform the fit for the two extremal values of the signal strength : /Uj| ex p[l ± 
<5/Xi//ij| t h] with the theoretical uncertainty <5/ii//Ui| t h conservatively assumed to be 20% for both 
the gluon and vector boson fusion mechanisms. 

The results of the fit of the data is presented in Fig. [I] and relies on established values [12] of 
Ax 2 = X 2 ~ Xmin with the following x 2 function 



/',('•/• 'V) - //, 12 



(J'i | exp 
i |exp 



In addition to the fit that is usually performed, leading to the colored regions in Fig. [T] for the 
best fits at the la, 2a and 3a levels in the c/ versus cy plane 8 , we also present the results of the 

fit for a treatment of the theoretical uncertainty as a bias. The value, /ij| exp + (5/ii//Xj| t h 

corresponds to the plain contours on Fig. [l] while the dashed contours are for the lower /i 
expectations (negative sign in front of the error). The distances between the two contours 
represent the theoretical uncertainty induced on the fitted parameters. 

A first conclusion that can be drawn from this figure is that, with the latest LHC data, in 
particular the new CMS di-photon rate which has no excess compared to (and is even slightly 
below) the SM expectation, the SM point is now included inside the la domain. Before the 
CMS update, there was a large deviation in the di-photon channel which triggered discussions 
about the possibility that c/ < (as shown in the figure for the 2a and 3a regions) which 
leads to a constructive interference between the top quark and H^-boson loop contributions 



[8] Because of the exact reflection symmetry under, c — > —c, leaving the squared amplitudes of the Higgs rates unaffected, only the 
part, cy > 0, is presented in Fig.^ 
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FIG. 1: Best-fit regions at 68.27%CL (green), 95.45%CL (yellow) and 99.73%CL (grey) in the plane c f versus c v , based 
on the % 2 function of eq. the best-fit location is indicated by a (black) cross. The 'concentric' best-fit domains at 
the same CLs obtained also from \ 2 but for the two extreme theoretical predictions - upper (plain contours) and lower 
(dashed contours) - of the Higgs signal strengths, are presented. The SM (red) point at c/ = cv = 1 is also shown. 



to increase the di-photon rate. This possibility has been elaborated within several effective 
scenarios in the recent literature to generate such a di-photon enhancement (see examples in 
Ref. [12]). 

The statement that the fit result is compatible with the SM expectation is true regardless of 
the method chosen to implement the theoretical uncertainty. Nevertheless, it appears clearly 
on the figure that treating this uncertainty as a bias enhances significantly its effects in the 
c-parameter determination. For instance, the entire la domain on the plot is larger when taking 
into account the possible shifts due to the theoretical uncertainty. 

As this theoretical uncertainty cancels out in the ratios Dxx of signal strengths of eqs. ( 4p ), 
we have performed a fit based on the Xr function 9 : 

Xr ^su 2 tew] 2 

We have considered the inclusive di-photon channels of CMS [U [29] and ATLAS [32] that 
are largely dominated the ggF mechanism. Regarding the ZZ final state in ATLAS [33] and 
CMS [38], we have also used inclusive production 10 . Finally, for the WW and rr searches 
in ATLAS [3H [36] and in CMS [391 E]> we have selected Higgs production in ggF with an 
associated 0/1 jet or the VBF production mechanism. Hence for both ATLAS and CMS, the 



[9] The fact that the distribution of fi ratios is not gaussian is not expected to modify our results. Also, we refrain from including 
the H — > bb channel as the ATLAS and CMS sensitivities are still too low l3l 1351 . 
[10] In fact, for fJ-yy/fJ.zz\c>Lp, we will use the more accurate p^^/zz value provided by the ATLAS collaboration [2]; this quantity 
corresponds exactly to the branching fraction ratio and it is deduced from a combination of the ggF and VBF channels. 
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FIG. 2: Left: Best-fit regions at 68.27%CL (green), 95.45%CL (yellow) and 99.73%CL (grey) in the plane c/ versus c v , 
based on the Xr function. The best-fit (dotted) contours obtained from the \ 2 function in case of a theoretical error 
added in quadrature (as in Fig. [TJ are superimposed (in red). The associated best-fit point (cross) and SM (red) point 
are also shown. Right: Same plot as the left one but the best-fit domains from the \ 2 analysis are now derived for the 
two extreme theoretical predictions of the signal strengths (as in Fig. [Tj . 



situation is equivalent in a good approximation to have vanishing selection efficiencies except, 
e 33 — e 99 — 1 an< ^ e gg — e ^ — 1 or e vBF — c vbf — 1, so that the theoretical predictions for the 
ratios simply read as in eq. (|5]). 
The combined ratio values measured by the ATLAS and CMS collaborations are 



fJ-ZZ 



_ 1 1 +0.4 ^TT 

- i - i -0.3 > 

cxp H>WW 



VBF 



-0.24 ±0.83 and 



exp ^WW 



99 



cxp 



1.2 ±0.75 (10) 



The errors 5(/^ 77 //izz) and 5(fi TT / \iww) are computed assuming no correlations between the 
different final state searches. These uncertainties on the ratios are derived from the individual 
errors, Sfii - provided in the experimental papers - and are thus also dominated by the exper- 
imental uncertainties, e.g. 8(fJL /y / \izz) ~ fK/Vr/ZAz^lexp, as expected from the fact that the 
theoretical uncertainty largely cancels in ratios _D 77 and D TT . These ratios are given by 



\\cwAi_\mw] + (|) 2 c f A[m t ] + {-\fc b A[m b ] + (§) 2 c c yl[m c ] + \c T A[r. 



\\Ax\mw] + (l) 2 A[m t ] + (-^A[m b ] + (|) 2 A[m c ] + \A[m T ]\ 2 

D TT ~ , (11) 

\cw\ 

where A[m] = 74 1 / 2 [r(m)] and Ai[r(m)] are respectively the form factors for spin 1/2 and 
spin 1 particles [8j normalized such that A\r{rn) <C 1] — > 1 and Ai[r(m) C 1] —7 with 
r(m) = Mjj/Am 2 and, for m H — 125 GeV, one has v4 1 [r(m l y)] ~ —8.3 and Ai/ 2 [r(mt)] w 1. 

In Figure [2j we show the results from fitting the Higgs decay ratios through the function 
Xr- In the left panel, the best-fit domains obtained e.g. at la do not exclude parts of the 
la regions obtained from x 2 ; such a compatibility was expected since the main theoretical 
uncertainty cancels out in the Dxx ratios and is negligible for the signal strengths since it is 
added in quadrature to the experimental error as already described. The domains from x 2 are 
even more restricted as (i) this function exploits the full experimental information on the Higgs 
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FIG. 3: Best-fit regions at 68.27%CL (green), 95.45%CL (yellow) and 99.73%CL (grey) in the plane c t versus c v as 
obtained from a three-dimensional fit (whose best-fit point is the black cross on the central plot) of the 3 free parameters, 
Cf, c t , cy, based on the x\ function. The three two-dimensional plots correspond to the slices of these three-dimensional 
domains at c/ = 0.52, 1.14 and 1.55. Superimposed are the best-fit domains at 68.27%CL, 95.45%CL, 99.73%CL obtained 
from x 2 for the two theoretical signal strength predictions - upper (red plain) and lower (red dashed contours). 



rates and not only on the ratios and (ii) the experimental error on a ratio of rates is obviously 
higher than on the rates alone. 

In the case where the theoretical error for each Higgs channel is taken into account as a bias, 
the best-fit contours span wider regions of the parameter space. This could be seen in Fig. [T] 
and the same contours appear in the left-hand side of Fig. [2] as well as in its right-hand side 
part where the theoretical uncertainty enters as a bias for \ 2 ■ I n the latter case, the large la 
domain from x\ excludes a small part of the la region (near the SM point) obtained from x 2 ■ 
Besides, the lower 95.45 %CL region from the x 2 fit is completely excluded by the x 2 r domain (in 
yellow) at the same CL. In conclusion, within the more realistic case of treating the theoretical 
uncertainty as a bias, the x\ domains can thus play an important role by excluding parts of 
the x 2 -fit regions; this is due to the increased contribution of the theoretical error, in x , which 
does not affect x 2 r- 

In Fig. [3j we present the three-dimensional x 2_ fit results in the case where the parameters Cf, 
c t and cy are free. It shows how precisely are presently known the top quark Yukawa and the 
gauge boson interactions with the Higgs scalar (in the case of a preserved custodial symmetry). 
For either lower (left plot) or higher (right plot) bottom quark Yukawa couplings as compared to 
the SM, the size of the characteristic la regions decrease significantly, which also gives an idea 
of the present knowledge of the coupling c&. This more realistic three-dimensional fit illustrates 
as well the potential interest of the X 2 r~^ '■ f° r instance, one observes on the central plot, that 
it excludes at la the lower-right (i.e. the dysfermiophilia solution c t < 0) and upper-left (i.e. 
its almost symmetric domain) la regions resulting from the X 2_ fit- 

Since x\ is only affected by the experimental uncertainty, it is interesting to quantify the 
evolution of the fit when the experimental systematic and statistical errors are reduced [43J. 
For this purpose, we combine the present measurements with the expected results from the 
14 TeV LHC in each channel investigated by ATLAS and CMS. We assume the central values 
at 14 TeV to be identical to those from the combination of the 7 and 8 TeV data, and that the 
future experimental errors, 5//j| exp , will reduce essentially like the inverse of the square roots of 
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FIG. 4: Best-fit regions at 68.27%CL (green), 95.45%CL (yellow) and 99.73%CL (grey) in the plane c/ versus c v , based 
on the Xr function and including hypothetical data from the 14 TeV LHC with C — 300 fb" 1 [left plot] or 3000 fb _1 
[right plot]. The best-fit A\ 2 contours (with their best central points) obtained in the same conditions for the two 
theoretical predictions of signal strengths are superimposed (red plain and dashed contours); the (red) bands joining the 
two la ellipses indicate typically the theoretical uncertainty (band length) and experimental error (band width). The 
exactly symmetric domains, obtained via c/ — > —Cf, c v — > —c v , are not shown. 



number of events, y/ojC with C the integrated luminosity 11 . 

The estimated x 2 fit results at 14 TeV are presented in Fig. [4] assuming luminosities of C = 300 
and 3000 fb _1 |43j. The behavior of the best-fit x\ regions appearing in the figure originates 
from the compensation between the enhancement of T(H— ^77) and that of F(H— >ZZ) as cy 
increases, leading to relatively stable _D 77 values; the increase of T(H^-tt) and T(H— >WW) 
with increasing c r or cw also compensate each other in D TT . The best-fit values of Cf and cy 
in Fig. [4] should be considered as an illustration since the precise central values are of course 
not known, neither the exact experimental uncertainties. However, the above estimation of the 
statistical error provides an indication of the typical relative sizes of the best-fit \ 2 and x\ 
domains in the future. 

The main features that the plots of Fig. [4] exhibit are that when increasing the luminosity and 
hence, reducing the experimental error as shown by the smaller ellipses on the right plot, the x 2 - 
fit reaches the level where the theoretical error is dominating and fixes the typical uncertainty 
scale (stable red band length on the two plots), whereas for the X 2 r~^ m which the theoretical 
uncertainty is absent, the precision obtained on the couplings Cf and cy improves as long as the 
experimental error decreases (decrease of the colored region widths on the right plot). Thus, 
for high LHC luminosities, the fit of the decay ratios will play a crucial role and will have to 
be combined with the common Higgs rate fit, as illustrated on the right plot of Fig. (H): there 
for instance the la region from x 2 (typically the red band) is wider than when restricted to 
its intersection with the x\ domain at la (green band). This corresponds, to an improvement 
of the whole accuracy from ~ 10% down to ~ 5% on both Cf and cy; with such accuracies 



[11] This is justified for the statistical error and corresponds to an optimistic situation for the systematic error, which is difficult to 
predict for each channel but depends partially on the background rate uncertainties which have a statistical behavior as well. 
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one starts to be really sensitive to deviations in the Higgs couplings arising in supersymmetric 
theories or composite Higgs models as, for instance, discussed in Ref. 



4. The parity or CP— composition of the Higgs boson 

As mentioned in the introduction, the observables such as correlations in Higgs decays into 
vector boson pairs [15J or in Higgs production with or through these states (16] that are usually 
used to probe the Higgs parity project out only the CP-even component of the HVV coupling 
even if the state has both CP-even and CP-odd components. Thus, in these CP studies, one 
is simply verifying, a posteriori, that a CP-even Higgs state has been indeed produced. The 
HVV coupling takes the general form (here, we assume cy > 0) 

<vv = c v (M v /v) gT (12) 

where cy measures the departure from the SM: cy = 1 for a pure CP-even state with SM-like 
couplings and cy = for a pure CP-odd state. Indeed, the coupling of a pseudoscalar A state 
to W/Z bosons is zero at tree- level and is generated only through loop corrections which are 
expected to tiny. The measurement of cy should allow to determine the CP composition of a 
Higgs boson if it is indeed a mixture of CP-even and CP-odd states. 

However, having cy ^ 1 does not automatically imply that the observed state has a pseu- 
doscalar component. As a matter of fact, the Higgs sector could be enlarged to contain other 
neutral Higgs particles Hi that have not been detected so far because they are too heavy or 
too weakly coupled. In this case, the sum of the squared couplings of each state H to gauge 
bosons, c Vi g 2 H wi should reduce to the SM Higgs coupling, g HV v Hence, cy < 1, could mean 
that there are other CP-even states which share the SM Higgs coupling to VV with the ob- 
served Higgs boson. Nevertheless, in all cases, the quantity fiyp — 1 — c v gi yes an upper bound 
on the CP-odd contribution to the HVV coupling 12 . 

In contrast to the couplings to massive gauge boson, the CP-even and CP-odd components 
of the state can couple to fermions with the same magnitude and one can write 



rrif 

9H ff = -i — 



Re (c/) + % Im (c f ) 75 



(13) 



where in the SM one has Re(cy) = 1 and Im(cv) = but in general, the normalisation of the 
coupling, Re(c/) 2 + Im(c/) 2 = |c/| 2 , should be taken arbitrary as in the previous section. 

Hence, one can consider the same effective Lagrangian as in eq. ^ where cy = cz = cw 
represents exclusively the CP-even component of the observed boson assumed to be one eigen- 
state of an enlarged Higgs sector. In contrast, the c/ = q, = c c = c t and c t parameters for light 



fermions and the top quark contain the CP compositions of eq. (13) with the possibility of a 
deviation of the normalisation |c/| 2 compared to the SM Yukawa interaction. 

In the case of the light fermions, one has Mh > mj so that chiral symmetry holds and the 
partial decay widths (the only place where they enter if one neglects their tiny contribution to 



[12] The best example of an extended Higgs sector with CP-violation is the minimal supersymmetric extensions of the Standard 
Model (MSSM) with complex soft-SUSY breaking parameters 1281 . One has then three neutral Higgs states Hi,H2 and H3 
with indefinite parity and their CP-even components will share the SM HVV coupling, c vl + c V2 + c v ^ = 1. There are no 
antisymmetric CP-odd couplings HiV^vV^" at tree- level and those generated at the one-loop level are extremely tiny |28| . 
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FIG. 5: Best-fit regions at 68.27%CL (green), 95.45%CL (yellow) and 99.73%CL (grey) in the plane l-c v versus Im(c t ) 
for |c t | 2 = |c/| 2 = l ; these regions are obtained from a two-dimensional fit based on the Xr function. Superimposed 
are the best-fit regions at 68.27%CL, 95.45%CL, 99.73%CL obtained from \ 2 f° r the two theoretical signal strength 
predictions (plain and dashed contours in red). The SM (red) point is represented at the origin together with the best-fit 
points (the two black crosses) for the \ 2 n ts- 



the loop induced vertices) can be simply written as T(H — > ff) oc Re(c/) 2 + Im(c/) 2 oc \cf\ 2 
and the discussion in the previous section should entirely hold. 

In the case of the top quark, the situation is different as m t > Mh- The top quark enters 
the Hgg and if 77 vertices and the loop form factors for the CP-even A^ 2 and CP-odd Af, 2 

parts are in principle different [8]. Fortunately, in these vertices the approximation m t 3> Mh 
is extremely good and in this limit, the form factors take simple forms: Af, 2 = | and Ay, 2 = 2. 

Ignoring the small contributions of the light fermions for simplicity, the Higgs rates normalized 
to the SM expectations can be written as, 

T(H -> 77 ) _ WcwA^mw] + (f ) 2 Re(c t )| 2 + | (|) 2 |lm(c t )l 2 
r(H^ 77 )|sM ~ ' \i Al [ mw ] + (%) 2 \ 2 

a(gg^H) _ F(H -> gg) „ + ( (M) 



a(gg^H)\ SM T(H -»■ gg)\ S M 1 1 2 

with ylifmiy] — —8.3 for Mh ~ 125 GeV. For a pure pseudoscalar state, Re(cf) = 0, there is no 
W contribution to the H — > 77 rate; there are also no H — y ZZ and WW decays, a possibility 
that is clearly excluded by the present data as the A£ and 2£2v signals have been observed. To 
quantify the degree of exclusion of this possibility, one needs to measure kqp = 1 — Cy ( an d 
ideally, independently of the fermion couplings c/). 

Based on these rates, one can perform a fit using the same \ 2 function as in eq. ^ but with 
the dependence, \ 2 — X 2 [^ e ( c t), I m (ct), c/, cy\. The numerical results are displayed in Fig. ([5J 
for the x 2- fit an d t ne Xirfit w hich reveals itself to be useful as well for measuring the CP-odd 
component of the Higgs boson. In this case, we have made the simplifying assumption (besides 
cy > 0) that the absolute normalisation of the fermion couplings is the same as in the SM, 
I Ct 1 2 = \cf\ 2 = 1, but in the case of the top quark, Im(c t ) is assumed to be free. We assume 
Re(ct) > and the obtained plot is symmetric under Im(ct) — > — Im(ct). 

The conclusion is that, at the 99.73%CL or at the 3cr level, the CP-odd component of the 
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FIG. 6: Domains at la (green), 2a (yellow) from the central value of \izz, in the plane cv (left plot) or Cf (right plot) 
versus rg v / r l M [r| M being the total SM Higgs width] for c t = 1 and c f = 1 (left plot) or c v = 1 (right plot). The 
dependence of these constraints on the theoretical uncertainties is illustrated by the (black) curves which indicate the 
other possible extreme domains. The direct upper limit on from invisible searches at LHC (derived for cy = 1) |26 
is shown on the right plot. 



observed Higgs boson obeys the upper bound kqp — 1 ~ c v < 0.68. A pure CP-odd Higgs 
state, i.e. the case kcp ~ 1; is excluded with more than 4a. This is much more severe than 
the constraint from the correlations in H — > ZZ decays, which (with its inherent limitation 
discussed above) allows only a < 3a discrimination between the CP-even and CP-odd cases 



5. The invisible Higgs decay width 

In the previous discussion, the signal strength fizz in the channel H — > ZZ — > 4£ ± played a 
prominent role because the theoretical ambiguities are minimised: the measurement is inclusive 
and does not involve the additional theoretical uncertainties that are introduced when breaking 
the cross section into jet categories and there is no loop induced new physics effect as in the 
H — > 77 case; besides that, it is the most accurate single signal strength measurement. One 
can also use \xzz for the determination of the invisible Higgs decay width which enters in the 
signal strength through the total decay width T^, fi Z z oc T(H — > ZZ)jY^ with 

r tot = r nv + r SM (c/)Cy) (15) 

r^ vl (cj,cy) is the SM total width which is calculated with free coefficients Cf and cy and 
including the state-of-the art radiative corrections [20, 23J. One can write the ZZ signal strength 
as a function of and the Higgs couplings, fizz\th = Hzz(XW ' i c v, Ct, Cf = c c = Cb = c T ), and 
impose that it lies within its la or 2a ranges. This restricts the parameter space to specific 
regions as is shown in Fig. [6] where in the left-hand hand side c/ is SM-like and cy is varied 
and in the right-hand side, it is the opposite c/ is varied while cy = 1. 

On the figure, we also display for comparison the recent direct limit [26] on the invisible Higgs 
width obtained from combining the 7+8 TeV LHC data in the qq — > ZH — > Z+E/p direct search 
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channel. This gives B™ < 65% at the 95.45%CL if the assumption cy = 1 is made. With the 
simplifying Cf — Cy — 1 assumption, the indirect limit on the invisible width that one obtains 
from the signal strengths is better, as can be seen from the figure; at la, it reads as 

yUinv = r^ v /r| M < 0.52 @68%CL (16) 

This limit is at least a factor of two worse than those obtained in the similar fits of Refs. [14J 
using the latest LHC data, the reason being that, here, we assume a 20% theoretical uncertainty 
on the Higgs production cross sections that we treat as a bias and do not combine quadratically 
with the experimental uncertainty. 

6. Conclusion 

We have analyzed the Higgs production cross sections at the LHC for the different Higgs 
decay channels that have been searched for, H — > WW, ZZ, tt, bb and H — > 77. Using the 
latest results given by the ATLAS and CMS collaborations with the ~ 25 fb _1 data collected in 
the runs at = 7 and 8 TeV, we have first performed a fit of the Higgs couplings to fermions 
and massive gauge bosons and shown that they are now compatible with the SM expectation at 
the la level. The accuracy of the various experimental measurements is now almost saturated 
by the theoretical uncertainties stemming from QCD. 

We have argued that ratio of cross sections times branching ratios in different Higgs search 
channels are essentially free from these uncertainties and do not require further theoretical 
assumptions, on the total Higgs decay width for instance. These ratios, in particular in the 
H — > 77 v.s. H — > ZZ and H — > tt v.s. H — > WW channels, are being measured quite 
accurately already with the present data and provide tests of the SM predictions in a less 
model-dependent way. We show that at the 14 TeV LHC with a high luminosity, 300 fb _1 
and even 3000 fb _1 , they could allow the measurement of ratios of Higgs couplings with an 
accuracy at the level of a few percent which should allow to test the small deviations expected 
in realistic new physics models. 

In a second part of this paper, we have considered together with the ratios of cross sections 
times branching ratios in the most important search channels, the signal strength in the ex- 
tremely clean H — > ZZ channel in which the theoretical uncertainty is taken to be a bias. 
We have then shown that first, the particle observed at the LHC is at most 68% CP-odd at 
the 99%CL and the possibility that it is a pure pseudoscalar state (and hence does not couple 
to ZZ states at tree-level) is excluded at the 4cr level when including both the experimental 
and theoretical uncertainties. The signal strengths in the H — > ZZ channel also measure the 
invisible Higgs decay width which is shown to be T^/rf^ 1 < 0.52 at the 68%CL if the Higgs 
couplings to fermions and gauge bosons are assumed to be SM-like. 

All these results give us great confidence that the state observed at the LHC in July 2012 
is indeed a Higgs particle and, more than that, it resembles very closely to the Higgs particle 
predicted in the Standard Model. 
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